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ABSTRACT 
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A metal-rich environment facilitates planet formation, making metal-rich 
stars the most favorable targets for surveys seeking to detect new exoplanets. 
Using this advantage to identify likely low-mass planet hosts, however, has been 

£3 \ difficult: until now, methods to determine M-dwarf metallicities required obser- 

vationally expensive data (such as parallaxes and high- resolution spectra), and 
were limited to a few bright cool stars. We have obtained moderate (R~2700) 
resolution K-band spectra of 17 M-dwarfs with metallicity estimates derived from 
their FGK companions. Analysis of these spectra, and inspection of theoretical 

^ I synthetic spectra, reveal that an M-dwarf 's metallicity can be inferred from the 

strength of its Na I doublet (2.206 /im & 2.209 /im) and Ca I triplet (2.261 
/im, 2.263 /im & 2.265 /im) absorption lines. We use these features, and a 
temperature-sensitive water index, to construct an empirical metallicity indicator 
applicable for M-dwarfs with near-solar metallicities (-0.5<[Fe/H]<+0.5). This 
indicator has an accuracy of ±0.15 dex, comparable to that of existing techniques 
for estimating M-dwarf metallicities, but is more observationally accessible, re- 
quiring only a moderate resolution K-band spectrum. Applying this method to 
8 known M-dwarf planet hosts, we estimate metallicities ([Fe/H]) in excess of 
the mean metallicity of M-dwarfs in the solar neighborhood, consistent with the 
metallicity distribution of FGK planet hosts. 
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Subject headings: stars: late-type — stars: abundances — stars: planetary sys- 
tems 



Introduction 



M-dwarfs are under intense scrutiny as potential planet hosts. M-dwarfs have sufficiently 
low masses (O.6M > M+ > 0.08M Q ) and small radii (O.6i? > i?* > 0.1R Q ) that exoplanets 
induce considerably larger reflex velocities and transit depths than an identical planet would 
around larger, more massive hosts. Accordingly, most of the lowest-mass planets detected to 
date orbit low-mass stars. M-dwarfs' relatively low luminosities also ensure that habitable 
planets will be in small, fast orbits, producing higher reflex radial velocity amplitudes , and a 



greater likelihood, depth, and frequency of transit events (INutzman fc Charbonneau 



M-dwarfs constitute more than two thirds of the nearby stellar population ( jChabrier 



2008). 



2003|); 



identifying the most likely planet hosts in advance would greatly enhance the efficiency of 
M-dwarf exoplanet surveys. 

Observations of FGK-dwarfs have convincingly shown that gas giant planet frequency 
rises steeply with host star metallicity, a result known as t he planet-metallicity correlation 
(JGonzaleall997l : ISantos et al.ll2004J : iFischer fc Valentil 120051 ). Determining if this correlation 
extends into the M-dwarf regime requires accurate metallicity estimates for ultra-cool stars. 
As M-dwarf spectra are dominated by chemically complex molecular features, a priori spec- 
tral synthesis does not well match the observed spectra. Most attempts to calibrate M-dwarf 
metallicity indicators have therefore focused on observational studies of wide binaries, where 
the robustly estimated [Fe/H] of the FGK-primary is assumed to describe the M-dwarf sec- 
ondary as well. This assumption will hold if both binary components formed from the same 
well-mixed molecul ar cloud, and no mass transfer or dredge-up has occurred in the system. 
Bonfils et al.l ( 120051 . B05 hereafter) adopted this technique to develop a photometric metallic- 
ity calibration that suggested tha t nearby M-dwarfs , including the planet hosts GJ 876 and 
GJ 436, were slightly metal poor. iBean et al.l (120061 ) also estimated subsolar metallicities for 
GJ 436, GJ 581 and GJ 876 from spectral synthesis of optical spectra. Though there were 
only a sma l l num ber of M-dwarf planet hosts in their samples, the results of the B05 and 
Bean et al.l (120061 ) studies were unexpected, since they imply that planet formation around 
M-dwarfs does not depend on the star's metal content. This finding is difficult to reconcile 
with the planet-metallicity correlation for FGK- dwarfs, and runs counte r to the predictions 
of the commonly accepted core accretion model (IPerri fe Cameron! 11974 ). 



Johnson fc Appd (120091 . JA09 hereafter) derived a new photometric metallicity calibra- 



tion using 6 M-dwarfs with wide, metal-rich FGK companions not present in the B05 study. 
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JA09 found that the previous photometric calibration systematically underestimated the 
metallicities of these metal-rich stars. JA09 attributed this discrepancy to the lack of metal- 
rich calibrators in the B05 sample and/or the use of visual magnitudes with poorly character- 
ized uncertainties. JA09 also estimated the metallicities of 7 M-dwarfs with planetary-mass 
companions, concluding that M-dwarf planet hosts are preferentially metal-rich, just like 
FGK planet hosts. JA09 suggested that the lack of Jovian planets around M-dwarfs is most 
likely due to the influence of stella r mass, not ste l lar m etallicity, on the planet formation 



process, consistent with the work of lLaughlin et al.l (J2004J ) on giant planet formation around 
red dwarfs. 



Schlaufman fc Laughlinl (120101 . SG10 hereafter) recently revised the B05 and JA09 pho- 
tometric calibrations, finding that they systematically underestimate or overestimate metal- 
licity at the extremes of their range. Their results hint that low-mass planets may be 
more likely to be found around metal-rich M-dwarfs, which disagrees with the almost equal 
probability of find ing Neptune-mass planets around a metal-rich or a metal-poor FGK-star 



flSousa et al.l 120081 ). 



Existing methods for estimating M-dwarf metallicities rely heavily on observations at 
visual wavelengths, where metallic absorption lines and molecules strongly influence a star's 
spectrum. Mid- to late- M-dwarfs are typically too faint for high-quality observations at 
visual wavelengths, wh ich has limited most meta llicity analyses to early- M-dwarfs and a 



few nearby stars (e.g. IWoolf fc Wallersteinll2005l ). Lacking a satisfactory way to estimate 



metallicities across the entire M-dwarf Sequence has limited our understanding of planet 
formation around the lowest mass stars. 

In this Letter, we present a technique for estimating the metallicity of field M-dwarfs 
from the strength of prominent K-band absorption features. We calibrated our technique 
using spectra of M-dwarfs with existing [Fe/H] estimates, adopted as a proxy for their overall 
metallicities, obtained from the analysis of a wide FGK-companion. Our technique does not 
depend on accurate parallaxes and V magnitudes, or high-resolution, high-signal-to-noise 
spectra, as required by previous techniques. Instead, our technique requires simply a K- 
band spectrum at moderate resolution, which can be efficiently obtained for many M-dwarfs 
with current near-infrared spectrographs. Our method will enable metallicity estimates 
to be generated for stars which are too cool or too distant for other methods, identify 
the most advantageous targets for planet-searches around cool stars, and thus improve our 
understanding of planet formation around low-mass stars. In §2, we describe our observations 
and data reduction. We present our analysis and results in §3, and summarize and discuss 
our conclusions in 54. 
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2. Observations and Data Reduction 

Near-infrared spectra of the stars in this work were ob tained with the TripleSpec spec- 
trograph on the Palomar Hale Telescope (IHerter et al.ll2008l ) during 7 observing runs between 
2007 and 2009. TripleSpec at Palomar has no moving parts and simultaneously acquires 5 
cross-dispersed orders covering 1.0-2.4 /im at a resolution of A/AA~2700. 

The spectra were reduced with an IDL-based data reduction pipeline developed by P. 
Muirheaqj. The data were sky-subtracted using a sky-frame made by median combining 
5 exposures with the science target placed at different positions along the slit. Each sky- 
subtracted exposure was then di vided by a normalized flat-field, wavelength calibrated and 
optimally extracted ()HorneJ|l986|) . The spectra were flux-calibrated and telluric corrected us- 
in g observation s of an A0 V star at a similar airmass with the IDL-based code xtellcor ^general 
by IVacca et~aD (J2003I ). 



Our sample consists of: 



Seventeen M-dwarfs with wide (> 5" separation), common-proper-motion solar-type 
companions to serve as meta llicity calibrators. The F GK-primaries have spectroscopic 
metallicity measurements bv lValenti fc Fischer! ([20051. S POCS Catalogue). The binary 
systems were s electe d from the Idiese fc Jahreissl (J1991I ) catalogue of nearby stars, the 
Poveda et al.l ( 119941 ) catalogue of ne arby wide binary a n d mu ltiple systems, and the 
list of new HIPPARCOS binaries bv JGould fc Chanamel (bo04h - 



Twenty-nine nearby M-dwarfs analyzed by JA09, to compare our spectroscopic results 
with their photometric [Fe/H]. 



Eight M-dwarf planet host s observabl e from Palomar Mountain: GJ 876 (jMarcy et al. 



1998 



Bonfils et al. 



Delfosseet al.l ll998f). GJ 436 (But 



er et al.l 12004!: ISantos et al. 2004), GJ 581 



GJ 649 J Johnson et al.l lioioh . HIP 57050 jHaehighipour et alJboiol ). and HIP 79431 



( lApps et al. 



2010|). 



200fih. GJ 849 (butler et alJl200fih. GJ 1214 Jcharbonneau et al1l2009h 



Twenty-nine M-dwarf spectroscopic standa rds defined by Kirkpatrick et al.l f 199lL KHM 
system) or with KHM spectral types from iHenry et al.l (119941 ) 



Figure [T] shows the {V-K^, M^-} plane for our sample. The JA09 calibration is only 
valid for stars with: 1) accurate Johnson V magnitudes and parallaxes from the Hipparcos 



1 http://www. astro. Cornell. edu/^muirhead/T^Downloads 
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Catalogue (Ivan Leeuwenll2008l ) , t he Yale Parallax Ca talogue (Ivan Altena et al.ll200lf ) or the 
TASS Mark IV Survey Catalogue JDroege et alJl2007h . and 2) V-K s colors between [3.9,6.6]. 
JA09 assume that the mean M# versus V-K5 relation in the solar neighborhood is equivalent 
to an [Fe/H]=-0.05 metallicity contour, with stars brighter than this being metal-rich, and 
fainter stars being metal-poor. Most of the mid-type and all of the late-type M-dwarfs fail 
both of these conditions: they lack accurate parallaxes and/or Johnson V magnitudes, so 
were not included in the JA09 calibration. 



3. Analysis and Results 

We have empirically calibrated several NIR spectral features to serve as effective M- 
dwarf metallicity indicators. Figure [2] shows the K-ba nd spectra of three M4- M4.5 dwarfs in 
our sample, along with updated PHOENIX models □ (IHauschildt et al.lll999l ) with effective 
temperatures equivalent to an ~M3.5-4 star and varying metallicities. The overall shape of 
the K-band spectrum is well represented by the updated PHOENIX models, and the Na I 
doublet (2.206 /mi & 2.209 /mi) and the Ca I triplet (2.261 /mi, 2.263 /mi & 2.265 /mi) are 
seen to strengthen considerably with metallicity. There is a discrepancy, however, between 
the strengths of these features in our spectra and in the models of comparable metallicity. 

We measured the equivalent width (EW) of the Na I and the Ca I features to differentiate 
between metal-rich and metal-poor M-dwarfs. Each spectrum was normalized and the EWs 
were calculated by computing the ratio of the area of a feature to a pseudocontinuum across 
that feature, computed from featureless regions on either side of it, using the the IDL-based 
function measure_ew by N. Konidaris and J. Harker. Uncertainties in the EWs were obtained 
from the errors in the estimate d pseudocontinuum and th e errors in the measured intensities 
of each feature, as described in lSembach fc Savage] (119921 ) . 



Figure [3^ compares the SPOCS [Fe/H] for our M-dwarf metallicity calibrators to the 
EWs of their Ca I triplets; Figure [3b compares the EWs of the Na I and the Ca I features 
for our M-dwarf sample. The correspondence between the Na I and Ca I features and the 
metallicity of M-dwarfs is evident: metal-poor M-dwarfs have lower EWs than metal-rich 
M-dwarfs. The strength of these features also depend on temp erature and surface gravity, 
howe ver, as do most atomic lines in the NIR spectra of M-dwarfs ( Ali et al.lll995t llvanov et al. 
2004J ). 



To explicitly account for the influence of stellar temperature on the strengths of Na and 



2 T. Barman, private communication 
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Ca lines, we include an independent, temperature sensitive spectral index in our analysis. 
K-ban d H^O absorption b ands been used extensively to diagnose M-dwarf spectral types 



McLean et al.l 120031 ); we di agnose the effective temperatures of stars in our sample 



(e.g., 

using the H 2 0-K index defined by ICovey et al.l (j2010l ): 



HoO-K 



(^(2.18 - 2.20))/(J r (2.27 - 2.29)) 
(^(2.27 - 2.29))/(J r (2.36 - 2.38)) 



where (J-"(a — b)) denotes the mean flux level in the wavelength range defined by a and b, in 
microns. The linear relation between H 2 0-K index and the KHM spectral types of spectral 
standards in our sample is: 



Sp.Type = 42.343 - 42.920(H 2 O-K) 
cx(Sp.Type) = 0.94 



(1) 



We report the spectral types estimated from Equation (1) for the metallicity calibration 
M-dwarfs and planet hosts in Table 1, rounded to the nearest half subtype. 

A star's surface gravity also influences its spectrum, but those are typically subtle 
effects. All o f our targets are dwarf stars, which reduces any potential gravity ambiguity in 
our analysis. iFernandez et al.l (120091 ) reported empirical values of log g~4.9 for five M-dwarfs 
in eclipsing binary systems, while values of log g~ 4.6-4.9 have been inferred from mass and 
radius values for M-dwarfs in iDemory et al.l ( 20091 ). 



To identify the best fit relationship between [Fe/H] , the EWs of Na I and Ca I features 
in A, and the H 2 0-K index, we performed a linear regression on the parameters measured 
for our M-dwarfs with SPOCS metallicities. We found a best-fit linear equation of: 



[Fe/H] = 0.142(NaI EW ) + 0.011(Ca I EW ) 
+ 2.541(H 2 0-K) - 3.132 
a([Fe/H]) = 0.15 



(2) 



with a residual mean square RMS( [Fe/H]) =0.02 and an adjusted squared multiple correla- 
tion coefficient i?^([Fe/H])=0.63, which implies, by comparison of the same criteria functions 
calculated in SG10, that our model provides a better fit than the models of B05 (RMS— 0.05, 
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Rl<OM), JA09 (i?M5=0.04, i?2=0.059) and SG10 (RMS =0.02, R 2 a =0A9). We estimate 
systematic errors, such as the zero-point value and the inclusion of a non-physically bound 
binary as a metallicity calibrator, to be of the order of ~0.06 (SPOCS data) and ~0.03, 
respectively. The combined systematic error is of the same order of the structure in the 
[Fe/H] residuals of the calibration sample. We note that the predicted [Fe/H] for some of 
the planet hosts lie above the [Fe/H] boundaries of the calibrators. The values of the EWs of 
Na I and Ca I, the H 2 0-K index and [Fe/H] values for the metallicity calibration M-dwarfs 
and planet hosts, are given in Table 1. 

The best fit linear combination of the EWs of the Ca I and the Na I features are plotted 
in Figure Hk versus the H2O-K index of our M-dwarf sample. There is a clear distinction 
between the M-dwarfs with metal-rich and metal-poor FGK-dwarf companions in Figure H^; 
the SPOCS [Fe/H] values of our M-dwarf metallicity calibrators are also plotted in Figure Hb 
against the value of their H2O-K index. The M-dwarfs with photometric metallicity estimates 
calculated by JA09 are also bisected between metal-rich and metal-poor using these K-band 
features. Moreover, we found that all of the M-dwarf planet hosts have metallicities much 
higher than the mean metallicity of M-dwarfs in the solar neighborhood ([Fe/H]sjv~-0.17), 
defined by SG10. The Jovian planet hosts are co-located with the M-dwarf companions of 
metal-rich FGK stars in the upper portion of Fig. Hk, while the Neptune planet hosts are 
located just above the M-dwarf companions of metal-poor FGK stars. 



4. Discussion 

We have found that the EWs of the Ca I triplet, the Na I doublet, and water absorp- 
tion in the K-band differentiate metal-rich and metal-poor M-dwarfs, as shown in Figure SJ 
We estimate [Fe/H] values higher than -0.05 dex for eight M-dwarf planet hosts, with the 
Jovian planet hosts being more metal-rich than their Neptune host analogs (Table 1). Our 
results, along with the conclusions of JA09 and SG10 that M-dwarf planet hosts appear to 
be systematically metal-rich, are consistent with the metallicity distribution of FGK-dwarfs 
with planets. 

We also estimated spectral types for the M-dwa rf planet host s and m etallicity calibrators 



in our sample, using the H 2 0-K index defined by ICovev et al.l (120101). O ur spectral type 
estimates agree well with the KHM spectral types ( JKirkpatrick et al.lll99ll ). modulo the ~1 
subtype errors (Table 1). We do find discrepancies, however, with the spectral types of some 
of the most metal-rich and metal-poor stars in our sample. The metal-rich planet hosts GJ 
849 and GJ 876 have later spectral types in the KHM system than are indicated by the H2O- 
K index. A possible explanation is that the KHM system relies strongly on the double-metal 



TiO and VO molecular bands in the 6300-9000 A wavelength range. Thes e molecular bands 
are known to be metallicity dependent (e.g. the subdwarf classification by lGizislll997l ): they 
grow stronger with increasing metallicity as well as decreasing temperature. While the KHM 
system includes the overall shape of the stellar spectrum, this does not change considerably 
for early M-dwarfs, such that metal-rich early-type M-dwarfs with strong molecular features 
could be easily mistaken for less-metal rich later-type M-dwarfs. This reasoning can also be 
applied to GJ 611 B, the most metal-poor star in our sample: we derive an M6 spectral type 
from its H 2 0-K index, while its KHM spectral type is M4, consistent with its weak TiO and 
VO bands being interpreted as a temperature effect rather than a sign of its low metallicity. 

Measuring M-dwarf metallicities will also improve our ability to determine key physical 
parameters for both stars and planets. The masses and radii of M-dwarf planet hosts are 
typically derived from theoretical models, and/or empirical relationships between these pa- 
rameters and t he star's temperature and luminosity, all of which have a stron g dependence 
on metallicity ( iBaraffe et al.l Il998l ; iDelfosse et al.l Il998t iDemory et al.ll2009l ). Since mini- 
mum planetary masses inferred from RV measurements are tied to the mass of the host star, 
and planetary radii derived from transit observations depend directly on the adopted stellar 
radius, errors i n stel lar parameters lead directly to errors in inferred planetary properties. 
Delfosse et al.l (119981 ) report a stellar mass of 0.3 M for GJ 876, but with an uncertainty 
of ~30% due primarily to the uncertainty in the host's metallicity, which affected the ini- 
tial estimate of the masses of its planetar y companions. Recent measured M-dwarf radii and 
masses have typical errors less than 10% (IFernandez et al.ll2009l ). but the errors may increase 
somewhat for stars with significantly non-solar metallicities. With poor estimates of the mass 
and radius of exoplanets, it is harder to constrain their true composition, atmospheres and 
interiors. 



Despite the fact tha t K-band alkali features are temperature/gravity dependent ( jAli et al 



1995c llvanov et al.ll2004l ). and prone to non-local thermodynamic equilibrium effects, our em- 
pirical analysis demonstrates that the dispersion in EW of the Ca I and Na I features can 
diagnose the metal abundance of M-dwarfs with roughly the same temperature (H2O-K in- 
dex) and gravityj. Our method has an advantage with respect to previous photometric and 
optical spectroscopic calibrations, since it does not depend on parallaxes or V magnitudes, 
allowing us to cover a larger sample of cooler and distant M-dwarfs. Our results enable 
the identification of the most advantageous M-dwarfs for planetary search surveys around 
cool stars, such as the Triplespec Exoplanet Discovery Instrument (jEdelstein et al.l 120071 . 



3 Future studies that may suffer contamination by giants should use other gravity sensitive features (e.g. 
the 2.296 /j,m CO band-head) to distinguish between dwarfs and giants prior to interpreting Na I and Ca I 
feature strengths as metallicity indicators. 
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TEDI), MEarth (jlrwin et al.ll2009[ ). the Pr ecision Radial Vel ocity Spectrograph (IJones et al. 



20081 . PVRS), and the CRIRES Survey bv lBean et all (120101 ). and will thus help unravel the 



formation mechanisms of planets around M-dwarfs. 
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HIP 79431 

/ GJ 849 
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GJ 876 




GJ 1214 



[Fe/H]=-0.05 



V-Kc 



Fig. 1.— The {V - Kg, M^} plane for our M-dwarf sample. The big red and blue dots 
are M-dwarfs in our metallicity calibration sample with [Fe/H]>-0.05 and [Fe/H]<-0.05, re- 
spectively. The black dotted line is a fifth-order polynomial fit to the mean Main Sequence 
calculated by JA09 ([Fe/H]=-0.05). The small red and blue dots represent M-dwarfs with 
[Fe/H] >-0.05 and [Fe/H] <-0.05 respectively, according to the photometric calibration by 
JA09. The big black dots represent the M-dwarf planet hosts. The lack of precise paral- 
laxes for most of the metallicity calibrators can explain the position of three metal-poor 
M-dwarf secondaries above the grey-dotted line, where metal-rich M-dwarfs are supposed to 
be according to JA09. 
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Fig. 2.- (a) K-band TripleSpec spectra of GJ 324 B (top), HIP 57050 (middle), and 
GJ 783.2 B (bottom), (b) Updated PHOENIX models with T e// =3200 K, log g=4.5, and 
varying overall metallicities. The regions us ed to calculate the EWs of the Na I and the Ca 
I features, and the H 2 0-K index defined by I Covey et al.l (120101 ) are shown in red, blue and 
grey, respectively The updated PHOENIX models are consistent with the overall shape of 
the M-dwarf spectra. However, there are discrepancies in the strength of absorption features 
between the M-dwarf spectra and the models of similar metallicity. 



-14- 



0.6 
0.4 
0.2 

O 

J? 

x 

aJ -0.2 

Ll_ 

-0.4 
-0.6 

-0.8 

10 

9 
8 
7 

ra 6 

-z. 

5 
4 
3 
2 



a) 



##H 






o 



© 



[Fe/H]=-0.05 



b) 



H h 



HtH 



** ^ 



H h 



HgH 



"* ^ 

-r T" _- ~ tAj — — — 



^^ ^H 







«#""■ • 



M4 



* h«H 



MO 



2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 



Ca I 



E\\ 



Fig. 3.— (a) SPOCS [Fe/H] versus the EWs of the Cal triplet for the M-dwarfs in the 
metallicity calibration sample. The nomenclature for the stars is the same as in Figure [TJ 
The dashed line is the mean metallicity of the solar neighborhood by JA09 ([Fe/H]=-0.05). 
(b) The EW of Na I versus the EW of Ca I features for our M-dwarf sample. The dashed 
lines are the [Fe/H]=-0.05 contours for the spectral types MO, M4 and M8, calculated from 
Equation (2). These alkali features correlate with [Fe/H], but the dispersion seen is partially 
due to differences in the stars' temperature: an M8-dwarf has larger EW values than an 
MO-dwarf of the same metallicity. 
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Fig. 4. — A linear combination of the EWs of the Ca I and Na I features versus the H 2 0-K 
index for our M-dwarf sample. The nomenclature for the stars is the same as in Figure [TJ 
The [Fe/H] values for the metallicity calibrators are also plotted versus the H2O-K index, 
to emphasize the index's insensitivity to metallicity. The dashed lines in the top panel 
are iso-metallicity contours for [Fe/H] values of -0.30, -0.05 and +0.20, calculated from 
Equation (2). According to our determination, the Jovian M-dwarf planet hosts have higher 
metallicities than M-dwarfs with Neptune- or Earth-size planets, which is in agreement with 
the metallicity distribution of FGK-dwarfs with planets. 



Table 1. M-Dwarf Planet Hosts and Metallicity Calibration Sample 





EW 


[A] 


Index 


This Work 


KHM 


SPOCS 


B05 


JA09 


SG10 


M sin(i) 




Name 


Nal 


Cal 


H20 


Sp. Type 


[Fe/H] 


Sp. Type 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


Mj 


Planet Notes 


HIP 79431 


9.699 


5.470 


0.904 


M3.5 


+0.60 


M3.0 






+0.16 


+0.52 


+0.35 


2.1 


Single Jupiter 


GJ 849 


8.043 


5.635 


0.955 


M1.5 


+0.49 


M3.5 






+0.14 


+0.58 


+0.41 


0.82 


Single Jupiter 


GJ 876 


8.126 


4.721 


0.930 


M2.5 


+0.43 


M4.0 






+0.03 


+0.37 


+0.23 


2.64 


2 Jupiters + Neptune + Super Earth 


GJ 1214 


8.520 


4.095 


0.895 


M4.0 


+0.39 


M4.5 








+0.03 


+0.28 


0.0179 


Single Super Earth 


GJ 649 


5.651 


4.722 


0.952 


Ml. 5 


+0.14 


M1.0 






-0.18 


+0.04 


-0.03 


0.328 


Single Neptune 


HIP 57050 


6.628 


4.410 


0.890 


M4.5 


+0.12 


M4.0 






-0.02 


+0.32 




0.298 


Single Neptune 


GJ 436 


5.328 


4.456 


0.915 


M3.0 


-0.00 


M2.5 






-0.03 


+0.25 


+0.10 


0.072 


Single Neptune 


GJ 581 


5.108 


4.202 


0.921 


M3.0 


-0.02 


M3.0 






-0.25 


-0.10 


-0.22 


0.0492 


Neptune + 3 Super Earths 


GJ 212 


5.067 


5.286 


0.963 


M1.0 


+0.09 


M0.5 


+0.19 


-0.05 




+0.09 










HD 46375 B 


6.830 


5.648 


0.934 


M2.0 


+0.27 


M1.0 


+0.24 


-0.33 


+0.22 


+0.12 










GJ 797 B 


4.654 


3.604 


0.933 


M2.5 


-0.06 


M2.5 


-0.09 


-0.09 




+0.11 










GJ 872 B 


4.621 


3.431 


0.930 


M2.5 


-0.08 


M1.0 


-0.22 


+0.12 














GJ 250 B 


4.893 


4.387 


0.925 


M2.5 


-0.04 


M2.0 


+0.14 


+0.12 




-0.07 










GJ 768.1 B 


5.677 


4.190 


0.921 


M3.0 


+0.06 


M3.5 


+0.16 


+0.16 




+0.22 








C 


NLTT 14186 


5.937 


5.564 


0.921 


M3.0 


+0.11 


M 


+0.05 
















HD 222582 B 


6.637 


3.919 


0.913 


M3.0 


+0.17 


M3.5 


-0.03 


-0.04 














GJ 777 B 


6.151 


3.177 


0.892 


M4.0 


+0.04 


M4.5 


+0.21 




+0.19 












GJ 231.1 B(C) 


5.439 


4.179 


0.892 


M4.0 


-0.05 


M3.5 


-0.04 


+0.16 




+0.12 










GJ 324 B 


7.867 


4.489 


0.892 


M4.0 


+0.30 


M4.0 


+0.31 




+0.33 


+0.26 










GJ 783.2 B 


4.716 


3.089 


0.883 


M4.5 


-0.19 


M4.0 


-0.15 


-0.35 




-0.19 










GJ 3348 B 


5.534 


4.154 


0.882 


M4.5 


-0.06 


M4.0 


-0.22 


+0.16 














GJ 544 B 


5.342 


3.605 


0.875 


M5.0 


-0.11 


M6.0 


-0.18 
















GJ 611 B 


3.105 


2.639 


0.856 


M6.0 


-0.49 


M4.0 


-0.69 
















NLTT 15867 


4.097 


2.905 


0.851 


M6.0 


-0.36 




-0.10 
















GJ 376 B 


8.418 


3.310 


0.821 


M7.5 


+0.18 


M5.0 


+0.20 
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